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Abstract 

Exclusive semileptonic B decays into excited charmed mesons (D\, D2) are 
studied up to the order of 1 /mg in the framework of the heavy quark effective 
field theory (HQEFT), which contains the contributions of both particles 
and antiparticles. Two wave functions rj^ and t]q, which characterize the 
contributions from the kinematic operator at the order of 1/ vtiq , are calculated 
by using QCD sum rule approach in HQEFT. Zero recoil values of other 
two wave functions and extracted from the excited charmed-meson 

masses. Possible effects from the spin-dependent transition wave functions 
which arise from the magnetic operators at the order of I/trq are analyzed. 
It is shown that the experimental measurements for the branching ratios of 
B — * D\ Iv and B — > D\ Iv can be understood in the framework of HQEFT. 
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I. INTRODUCTION 



Recently, studies on the semileptonic B decays into excited charmed mesons become 
interesting in both experimental and theoretical sides. Experimentally, to precisely measure 
the branching ratios of the semileptonic B meson into the groundstate charmed mesons, 
it also needs to measure precisely the branching ratios of the semileptonic B decays into 
excited charmed mesons which is the main backgroud for the former decays. Theoretically, 
it provides additional modes for testing the validity of effective theories, in particular, how 
good of the spin-flavor symmetry. 

The semileptonic B decays into excited charmed mesons have been discussed in [|IJ.[2|] 
based on the framework of the usual heavy quark effective theory (HQET). Where the de- 
pendence of decay rates on wave functions was presented in a general form, the authors 
mainly deal with Isgur-Wise type function and the 1/rriQ order wave functions T\ and r 2 
which arise from the effective currents. The wave functions rff arising from the chromomag- 
netic term in Lagrangian have been neglected, and those arising from the kinematic term in 
Lagrangian have been absorbed into the Isgur-Wise type function and not been considered 
separately. The detailed field theoretic calculations on these decays were presented in [[3],|J, 
where the Isgur-Wise type function and r l5 r 2 have been evaluated by using the QCD sum 
rule approach. 

The purpose of this paper is to study the transitions between ground and excited state 
heavy-light mesons within the framework of heavy quark effective field theory (HQEFT) 
in which both quark and antiquark effective fields with keeping quark- ant iquark coupled 
terms have been considered. In particular, we focus on the semileptonic B decays into the 
jf = | + charmed meson doublet (Di, D^), where ji and P are the spin and parity of the 
light degrees of freedom in the charmed mesons. It has been shown that the HQEFT can 
provide a consistent description on both exclusive and inclusive ||,[| decays of heavy 
hadrons. It has been seen in these references that the HQEFT has many advantages with 
respect to the usual HQET. In the new framework, the values of \V c t,\ extracted from both 
exclusive and inclusive decays show a good agreement; the bottom hadron life time ratios 
can be well understood; 1 / rriQ order corrections at zero recoil in both exclusive and inclusive 
decays automatically absent without imposing the equation of motion iv ■ DQ V = when 
the physical observables are presented in terms of heavy hadron masses; less wave functions 
are needed; and there are interesting relations among meson masses and transition wave 
functions at zero recoil. In the most recent paper 0], the decay constants and binding 
energy of ground state heavy-light mesons as well as transition wave functions among them 
were studied consistently by using QCD sum rule approach in the framework of HQEFT. It 
was noticed that the HQEFT appears much more reliable than the usual HQET in describing 
the decays and transitions for the ground state mesons. Particularly, the 1/rriQ corrections 
to the heavy-light meson decay constants were found to be much smaller than the heavy 
quark mass, so that the scaling law of the decay constants is only slightly broken. This 
observation is unlike the usual HQET, which may lead to either the complete break down of 
the 1/rriQ expansion |L0| or a large 1/rriQ correction |Tl[] in evaluating the decay constants. 
Two transition wave functions Ki(y) and ^2(2/) at the l/mg order have also been evaluated 
in |7]]. Their zero recoil values agree with those extracted from the ground state meson 
masses. 
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Our present paper is organized as follows. In Sec.|I|, the weak transition matrix elements 
relevant to B — > (D^D^) decays are studied up to the order of 1/m.Q and parameterized 
by independent universal wave functions. The form factors and decay rates are then given 
in terms of those wave functions. The relevant normalization of the excited mesons is also 
discussed. In Sec. [IT], using the appropriate interpolating currents for excited heavy mesons, 
we derive the sum rules for two of the important wave functions concerned at I/tuq order, 
t]q and t]q. In Sec.[TV], we present our numerical results obtained from the sum rule approach. 
The B — > (DijDz) decay rates and branching ratios are discussed in detail. Finally, we 
come to our brief summary in Sec.|V[ 



II. ANALYTIC FORMULAE FOR B — ► {D\,D%) IN HQEFT 

The matrix elements relevant to the semileptonic decays B — ► (£)]_, Dg) are the ones of 
vector and axial vector currents {V^ 1 = 07^6 and N 1 = C7 M 7 5 6) between B and the excited 
doublet [Di,D^). Usually, these matrix elements are parameterized as 

(D l {v' 1 e)\V^\B{v)) 
{D 1 {v' ,e)\A^\B{v)) 
(D* 2 (v',e)\A»\B(v)) 
(D*(v',e)\V»\B(v)) 

where the form factors /, and fcj are dimensionless functions of y — v • v', and e* M (e* Ma ) is 
the polarization vector (tensor) of D\ (D^)- 

In the framework of HQEFT (its main formulation is presented in Appendix |A], we 
can in general introduce an effective heavy hadron state \H V ) for the witness of exhibiting 
manifestly the spin-flavor symmetry |6|-p. It is related to the hadron state \H) in the full 
theory via 

- 1 (H'\Q'TQ\H) = 1 _ (H'JJeffe'S^-"^). (2.2) 
y/™H>m H y/A H >A H 

with \H V ) being normalized as 

{H v \Q+ 1 »Qt\H v ) = 2kv» (2.3) 

where 

A = Atf - 0{l/m Q ) = lim A H (2.4) 

is taken to be heavy flavor independent, and it mainly reflects the effects of the light degrees 
of freedom in the heavy hadron characterizes the off-mass shellness of the heavy quark 
within the heavy hadron. Once A is chosen to be the flavor independent binding energy of 
the hadron \H), one yields 

A H = m H -m Q . (2.5) 



= y/mEjnElfv^ + (jVX + /y 3 ^)(e* • v)], 
= i yJm Dl m B fAt^ 1 e* v pv ^ , 

= ^m D *m B [k Al e*^v a + (k A2 v» + k A3 v'^)e* al3 v a v% 

= iJm D *m B k v e^eiyv p v' (2.1) 
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which is the total binding energy of hadron. 

The hadronic matrix elements can be expanded according to the order of l/niq. By 
including the I/tuq order corrections which can arise from both the current expansion and 
the insertion of the effective Lagrangian, one gets 



^L {0 \q T Q\M) -> (0\gTQ+\M v ) - -^(0\qT±.(i^Q+\M v ) + 0(l/mj) , 



A w ' Am <m'|qtq|m) -> (M^trg+K) - ^(M;,|Q:tr4r(^±)U + W 



^ m M tm M 2m Q iip\\ 

-^-«I»H ^) 2 ^rQ+|M v ) + 0(l/m 2 Q ). (2.6) 
2m Q' -i ^|| 

The form factors /$ and fcj can be parameterized by a set of universal wave functions. It 
is simplest to do this by using the trace formulism [0,EJ. The spin wave functions for the 
jf = | ground state mesons B, B* and jf = | + charmed mesons D\, D\ are 

. . _ /=■ f — 7 5 , for pseudoscalar meson , . 

M» — VAr + < , „ (2.7) 
I for vector meson 



^ = { - V% V ^ " 57^(y - «")], for Di (2 g) 

+ \ ^Ju, for D| 

The matrices .M„ and .7-^ satisfy the properties j>M. v = M. v = —A4 V ^, i)T% = T% = —T^i 
and J^7 M = T^v^ = 0. 

With respect to eq. ( |2.3|) and the normalization of the | + excited states, we parameterize 
the matrix elements in eq.( [2.6|) as follows 

(H v \Qt-fQ+\H v ) = e(y)Tr[^T va ], 
(H v \Q+^-^DlQt\H v ) = -^(^i-Trf^V^], 

{H v \Q^-^ l - aap F a PQ+\H v ) = ^Tr{A^ a \v,v')T va ^P + l -a a(i T va \ (2.9) 

with 

A aa ' a P(v,v r ) = a x {g aa g a>fi - g afi g a ' a ) + ia 2 g aa ' a afi + a 3 g aa> \v a ^ - v^ a ) 

+ ia A {g (Ta a a ' p - g afi (j a ' a ) + a 5 (g™v - g aP v a )Y' + ia 6 {g <T ' a a CTl3 - g a ' l3 cr CTa ) 
+ a 7 (# CT 'V - /V) T CT + a 8 a CT '<V a/3 + w 9 / ff ( 7 a / - 7 V). (2.10) 

Finishing the trace calculations, we get from ( |2.2| ) and (|2.3| ) 

2m Dl v» = ^{-2A'C + — «(1) + 54(l))}(e* • e)^, 
A Dl m c 

2m D *y = ^{2A'e' - — K(l) - 3/4(l))}e* ra e <ra i;'*, (2.11) 
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where 



/"» ■-) 



■(di — a 2 — CL4 — as — a$) 



(2.12) 



Eq.(gTTD yields 



A Dl = A' («i(l) +54(1)), 



711. 



Ad* = A' «(l)-3/c' 2 (l)). 

2 m r 



(2.13) 



In deriving these formulae we have used the normalization of the leading order wave function 



e'(i) = 1, 



(2.14) 



which is a direct result of eqs.fl2.4p and fl2.11| ). 

Eq. (|2.13f) is similar to those relations for ground state mesons @,^]: 



Ad = A («i(l) + 3K2(l)), 

m c 

Ad*=A-— (k!(1)-k 2 (1)). 



(2.15) 



For B — > (Z>i, D|) decays, we parameterize the relevant matrix elements as follows 



(H v ,\Q'jrQ v 

(H v >mr-^±-DiQt 

iv ■ D 
-iv'- D 

1 -iv'- D 



'+rn+\ Hv) = T (y)Tr[v j*,TM v ], 

1 



H v ) = -Vo^jTriv^TMv], 
H v ) = -v c (y)^Tr[v a ^,TM v }, 

H v ) = -iTriR^iv^^TP+ia^M, 
H v ) = -LTr[R^{v )V ')F°,ia a Pp' + TM 



(2.16) 



where R aa n and R^ap can generally be written as 

Kal3( V , V ') = vWlalp + V2 v aV' a lp + VsdvaVp, 
K a p( v i v> ) = VlV*lalf3 + V2 v v v a1p + VaQaaVp- 

t, rj\ and r/i are Lorentz scalar wave functions of y. 

Now finishing the trace calculations, and taking into account eqs. fl2.2Q , ( |2.1|) , (|2.6j ), (|2.13|) 
and ( |2.15| ), one finds that fi and ki are related to the transition wave functions as follows 



(2.17) 
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where 



= (1 " y 2 ){\r + 2^a(«i(1) + ^2{l)) + + 5 ^(1))] 



1 b 1 r - 3 



VS/v* = - 3 [f + ^(^(l) + 3«a(l)) + ^T7(«i(l) + 5^(1))] 



2mftA 2m c A 



= (y - 2)[f + 2^a(«i(1) + 3«2(1)) + ^jM'^ 1 ) + ^(i))] 

+ ~{y - 2)v b + + v) - - y) - % c (i - f y)], 

VQf A = -(1 + y){[f + 2^aM1) + 3«a(l)) + + 5 4(1))] 

fc^ = "(1 + y){[f + 2^x(«i(l) + 3«a(l)) + ^X7K(!) " 3 4(1))] 



= [f + 2^A (Kl(1) + 3K2(1)) + 2^ (K ' l(1) " 3K ' 2(1))] 

k v = - [f + ^Mi) + 3. 2 (i)) + ^K(i) - 34(D)] 



mfcA m c A' 



r = r 



2m b A 2m c A'' 

»7 6 = -3»7i - (1 - S/)»72 - \vi (2-19) 



III. QCD SUM RULE EVALUATION FOR r, rfi AND rj§ 

In order to calculate the wave functions r and t/q, t/q, one may study the analytic prop- 
erties of the three-point correlation functions 

i 2 J ^ze^*"^ (3.1) 

i 2 /^*e^^ (3.2) 
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* 2 j d 4 ^^'-*-^ 

= E*\u,u',y)L^ A , (3.3) 

?j*x*zeW ■ x - p - z) (0\T{J^ +t3/2 (x), (QT^DlQ)(0),4 t _ tl/2 (z)}\0) 

= 3*{u,u/,y)L%$, (3.4) 

* 2 / dVze^'^(0\T{.r w/2 (x),(Q D\ — Q^rQ)(0), J f _ 1/2 (*)}|0> 
7 -it)'- £) 

= S«(a;, t y > y)L^ 4 , (3.5) 

* 2 / dWze^' ^)(0|T{J 2 ^ 3/2 (a;),(Q S* — ^rQ)(0), J L >1/2 (z)}|0) 
J -it)'- D 

= Efi(u>,e/,y)I%$ (3.6) 

with T being 7^ and 7^7 5 for vector and axial vector heavy quark currents separately. 
Ly > jf wt ® are Lorentz structures associated with the vector and axial vector currents (see 
Appendix [B]) . 

Jj^Pji are the proper interpolating currents for the heavy-light mesons. Here 

4,- 1/2 = \J~~^Ql 5( l f° r pseudoscalar meson, 

j\ a _ jy 2 = ]j^Q'J±q for vector meson, (3.7) 

for I ground state doublet, and 

4+,3/2 = ^fQM-OPi - for I>i, 

4^3/2 = ^Qt^dlDi + llDI - 2 -gfiPi_)q for 23*. (3.8) 

for I doublet. 

Generally, the proper current Jj t p,j t for the state with quantum numbers j, P, ji have 
been investigated in [13, 3[]. These currents were proved to satisfy the following conditions 

(M 4,7-' Wf. I >'..,;) = ifp jl 6 jj ,6p P .6 jlj( rf*™ a > t 
i(0\T(J^(x) J^;J'(0))|0) = 5 jf S Ppf 5 jlfi (-iySgT Pl • • 

x J dt5(x - vt)U PJl (x) (3.9) 

in the limit mg — ► 00. r/" 1 ' " Qj is the polarization tensor for the spin j state, g°^ = g al3 — v' a v'^ 
is the transverse metric tensor, and 7" = 7° — v a (v ■ 7). 5 denotes symmetrizing the indices 
and subtracting the trace terms separately in the sets (cti • ■ -<x,) and (/?i ■ • •/%), /p,^ and 
IIpjj are a constant and a function of x respectively, they depend only on P and j\. 
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Eqs.( |3.9|) implies that the sum rules in HQET (or HQEFT) for decay amplitudes derived 
from correlators containing such currents receive contributions only from one of the two 
states with the same spin-parity (j, P) in the mg — > oo limit. And starting from the leading 
order, the 1/tuq corrections to the decay amplitudes can then be calculated unambiguously 
order by order. 

It should be noticed that the HQEFT differs from the HQET only in the 1/toq correc- 
tions. Therefore fpj t are the same in the two frameworks. For the ground state | mesons, 
the sum rule for /_,i/2 is also known ||10| , [14| in the HQET. It was also checked again in the 
HQEFT (where F = v^/-, 1/2) §] . Our result is 

A 1/2 e" 2A/T = ^ r duie-'^ - M(l + 4 ^ 



167r 2 7o 2 3tt 

(go^Fgg) 4a (a s FF) 



(3.10) 

For the | T doublet, the sum rule for /+,3/2 is found to be |13[] 

f +m e =M?J Q d ^ <" + * 12 32^ = ^ + ' 3/2 - (3 - U) 

The leading order wave function r for B — > {D\,Dl) decays is evaluated in |§ through 
studying the three point correlation functions ( |3.1| ), and ( |3.2| ). The Borel transformed sum 
rule for it reads 

+ l 12T 967T (y + iy~ Snr - (3 ' 12) 

The total external momenta in (|3.1| )- (|3~6| ) are p = vtiqv + k and p' = rriQ'v' + k' with k 
and k' being the residual momenta of the heavy quarks. 5' ? o(a;, uj', y) and 'E' n o(uj,uj',y) are 
analytic functions of uj = 2v-k + 0(l/rriQ) and uj' = 2v' -k' +0{l/m'o) with discontinuities for 
their positive values. Saturating the correlators in eqs.( [3.3|) -( [3.6|) with physical intermediate 
states in HQEFT, phenomenologically we represent them as follows 

(o|^t/ 2 l M/ )( M 1Q r ^^i^l M )( M l J o t ,-,i/ 2 l°) 

phen^VyA ~ (2A - a; - ie) (2A' - uj> - ie) 

+ / di/diA /V»/>X) -L£*« + subtractions, (3.13) 

Jd (z/ — u; — «e) (z/ — uj ' — ie) 

^ T ^ aP) v (QI^I^X^KQ gi I ^rg )(0) |M)(M| jt., /2 |o) 

V ' A m\m (2A-uj-ie)(2A'-uj'-ie) 

+ / dvdv'-, Pphys }^ V>) -L^/^ + subtractions (3.14) 

Jd {v - uj - ie){p' - uj' - ie) V ' A v ; 

with the first term in each equation being a double-pole contribution, and the second repre- 
senting the higher resonance contributions in the form of a double dispersion integral over 
physical intermediate states in the proper integration domain D. 



^% T fia(fiaf3) _ 
^ v hen Lj V,A — 



s 



The matrix elements in Q3.13| ) and (|3.14j) can then be transformed into series of the 



matrix elements in the HQEFT in powers of 1/tuq through (2.6). With (|3.9|) , the first 
terms in ( |3.13| ) and ( 3.14j) become in the limit rriq — > oo 



^rj b j-y.ainaP) _ /+,3/2/-,l/2 ~Vq r ^(jtap) /q -, p-\ 

-voieW,A ~ (2A-u>-ie)(2A'-u'-ie)~ Lv > A [ j 



^Vo T fia(^ial3) _ /+,3/2/-,l/2 r fMx(na/3) /q , fi x 

" (2A-o;-ie)(2A'-a/-ze)T^ A 



- / dvdv' ,^ pert — + E NP + subtractions (3-17) 



On the other hand, the correlation functions may be written as 

b(c) 

_„„„ = I dvdv'- , 

J-oo [y — uj — ie)(y' — uj' — ie) 

with the first term being perturbative contributions and the second non-perturbative ones. 
These can be calculated order by order in the framework of HQEFT by using the pertur- 
bation theory and the operator product expansion (OPE) as well. The Lorentz structure 
L^*|' ia ' 3 ' is extracted out and not presented in ( 13.171) . 

Assuming the quark-hadron duality, the sum rules for r/g and r/g read 

ri b(c) [ Ul [ V2 i P^ert 

H° ; _ = / dp I dp-. tt^ ; - + En p + subtractions. (3.18) 

p Jo Jo (p — oo — ie)(p' - uj' - ie) 

The next step of sum rule method is to perform the Borel operator 

BP = T lim -^t(-— )™ with T = — fixed (3.19) 

n^oo,-u;-»oo T(n) duj U 

to both sides of sum rules. Since there are two momentum variables uj and uj' for the 
correlators (|3.3| )- (|3~Jj ), a double Borel transformation B^f 'B^ should be performed to them. 



This has the effect to suppress higher resonance contributions on one hand, and to enhance 
the importance of low dimension condensates on the other, and thirdly, it also eliminates 
the subtraction terms. 

In QCD sum rule analysis for B semileptonic decays into ground state charmed mesons, it 
was argued [|ilj-|T6|| that the hadronic and perturbative spectral densities can not be locally 



dual to each other. But if one integrates the spectral densities over the "off-diagonal" 
variable z/_ = keeping the "diagonal" variable v + = fixed, the quark-hadron 

duality is restored in v + for the integrated spectral densities. This method was also used in 
to calculate the transitions between ground state mesons in the HQEFT. 

In the present case, the initial and final states belong to different doublets and are 
asymmetric. If one uses an asymmetric triangle in the perturbative integral, however, the 
resulting wave functions or their derivatives will unfortunately be divergent at y — 1 0. 
For these reasons, here we shall follow the method used in namely taking t' = t = 2T 
and integrating first the spectral density over v_ in the region —p + (p^(p + , we then obtain 
from the quark-hadron duality the following form 



*(i) 

riu, 

-'pole 



J Q dp + e- u+/T p pert (u + ) + ~ NP , (3.20) 
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where H denotes the result obtained by applying double Borel operators to E, and 

Ppert{v+) = dv^p pert (v + ,V-). (3.21) 

J — v+ 

In the present calculations, we will neglect the light quark mass and higher radiative 
corrections. For non-perturbatve terms, we consider only the contributions of the quark 
condensate, the gluon condensate and the mixed quark-gluon condensate. The relevant 
Feynman diagrams are presented in Fig.l. From the sum rules, we arrive at the following 
results 

, - Vo _(A+A')/T_ 1 1+% [»c W -q, + /T ,4 

a s (FF) 7-y 2a s (qq) T 2 

96vr {1 + y) 3 3tt {1 + y) 2 b 

(i) 



+ 3/2 -1/2TC 1 ^ ;/ = -- — / du + e +l UJl 



1 + 

Qs (FF) 9y + l _ 2^(3 + 2^ = 
96tt (1 + y) 3 3tt (1 + 2/) 2 V ' 

In the numerical calculations, we take the following typical values (a s = g 2 /47r) for the 
condensates 



(qq) « -(0.23 GeV) 3 ; 

i(qa al3 F a/3 q) « -m 2 (qq) with m 2 = 0.8GeV 2 ; 

a 

From Eqs.(|3~T0| ) -(|3~T2| ) , (|3722|) , one easily gets 



a s (FF) = ot s {F« Ff) « 0.04 GeV 4 . (3.23) 



r = - 




57e T 






-,1/2 x <S7^ 


-,3/2 






sn b 




A 




-,i/2 x sn^ 


-,3/2 






sn c 




A' 




.,1/2 x 5^ 


-,3/2 



(3.24) 

The QCD higher order radiative corrections have not been included in eqs.( |3.10"D -( |3.12D 
and (|3.22j ). But what we will use in our numerical analysis is eqs. (|3.24|) , which are ratios of the 
three-point correlators to the two-point correlators. Though the QCD radiative corrections 
may be large, one may expect that they may not influence the ratios in eqs.( |3.24"D significantly 
due to the cancelation of numerators and denomenators. 



IV. NUMERICAL ANALYSIS OF WAVE FUNCTIONS AND DECAY RATES 

Now we turn to the numerical analysis of the sum rules obtained in the previous sec- 
tion. Imposing usual criterium that both higher order power corrections and the contribu- 
tions of the contimuun should not be very large, we find the proper sum rule "windows": 
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0.7GeV(T(1.2GeV. f-,1/2 an d f+,3/2 have been studied in [jTWTOlfl and the corresponding 
thresholds were found to be in the ranges 1.6GeV(o; c (2.2GeV and 2.7GeV(cj**(3.2GeV. In 
Fig. 2 we present r(l), T]q(1) and ?7g(l) as functions of T at fixed values of cu* and uj*^ 1 ', 
where uj c = 1.9GeV and u** = 2.95GeV are used. We find that r(l) are stable around the 
threshold u* ~ 2.35GeV, whereas r/o(l) and t)q(1) are stable around a smaller threshold value 
uj*W pa 1.85GeV. With this analysis, we obtain the following values for the wave functions 

r(l) =0.8±0.1GeV 



A 



A' 



0.35 ± 0.04GeV 
1.15±0.15GeV (4.1) 



Where the errors mainly arise from the threshold. In Fig.3, the variations of r, i]q and t?q 
with respect to y are shown, where we have used T = lGeV. 

We now come to consider the B — ► (Z?i,Z)|) decay rates. The differential decay rates 
are given by 

dT(B - D 1 W) GUV^ml^ , fnH n __ , 2Uf3 , 2 iw2l 



<&/ 487T 



3 r l 



y2-l{2(l-2yr 1 + rf)[^ i + (^-l)/i] 



+ [(y-ri)/ Vi + (y a -l)(/ % +riA & )] 2 } l 

W - 1) 3 / 2 {3(1 - 2yr 2 + r 2 )^ + (y 2 - 1)A; 2 



1447T 3 

+ 2[(y - r 2 )fc Al + (y 2 - 1)(A: A 3 + r 2 k M )] 2 } (4.2) 

with the kinematically allowed ranges l(y(1.32 for B — > DiW and 1(^(1.31 for B — > D* 2 lv. 
The form factors /j, fcj are related to the wave functions as shown in (|2.18| ). The zero recoil 
values of K\ and k 2 in ( |2.18 ) have been evaluated by fitting from the ground state meson 



masses and also by the QCD sum rule method 0. Similarly, «^(1) can also be extracted 
from fitting the meson masses given in (|2.13 ), i.e., 



171 

4(1) = — -(mm -m Dl ). (4.3) 



To extrat we consider the spin average mass of each doublet: 



with 



n^m H _ + n + m H+ 

m H = ; (4.4) 

n_ + n + 



1 

n± = 2j± + 1, j± = ji ± -. 



( [H3I ), ( FTSD and O yield 

[(m' B -m B ) - (fh' D - m D )]m b m c 



«i(l)-«!(l) 



m 6 - m c 
m b -m c 
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Here rriB(D) are the average masses of the jf = | ground state mesons, whereas rn' B , D s are 

the ones of the jf = § + doublet mesons. 

When taking the average meson masses to be rhp = 1.971GeV, m' D = 2.445GeV, % = 
5.314GeV and m' B = 5.73GeV and the quark masses mj, = 4.8GeV, m c = 1.35GeV, we 
arrive at the following values 

« x (l) = -0.56GeV 2 
k 2 (\) = 0.05GeV 2 
«i(l) = -0.83GeV 2 

4(1) = 0.00675GeV 2 (4.6) 

rjf characterize the matrix elements of the chromomagnetic operator. They are often 
neglected from the argument that the chromomagnetic operator must have small effects due 
to the small D\ — D\ mass splitting M. In this section, we first neglect them but discuss 
their possible sizable effects late on. When rjf are neglected, the formulae of decay rates 
become very simple because each differential decay rate depends only on a composite wave 
function. 

dT{B ^ y DlW) = ^^ rliy + l)(y 2 - lf\y - 1)(1 + nf + y(l - 2 yn + r 2 )]f 2 ( y ), 

dT(B -> D*W) _ G 2 F \V cb \ 2 m% 
dy ~ 72^3 

with 



4(y + 1)(!/ 2 - Iplfe + - r 2 ) 2 + 3y(l - 2yr 2 + r 2 2 )]f 2 2 (y) (4.7) 



,b 



It is seen that the 1/mq corrections from the kinematic term do not change the relative values 
for the two differential decay rates. As the spin-symmetry breaking term k' 2 arising from the 
mass splitting is small, without including the l/rrig corrections from the chromomagnetic 
terms, the relative value of their total decay rates should have the similar behavior as the 
one at the leading order, i.e., 

T^B -> D x lv) ~ r x (5 -> D|/P) 1 J 

which can be explicitly seen from Table 1, where we have used the input values rrib = 
4.8GeV, m c = 1.35GeV, \V c b\ = 0.038, the life time tb = 1.6ps and the thresholds u c = 
1.9GeV, u* c * = 2.95GeV, to* = 2.35GeV, to*^ = 1.85GeV. Note that both decay rates of 
B — ► DiW and B — > receive large I/tuq contributions (but not as large as the 1/tuq 
contributions received by the B — > Di/z/ decay rate given in |17J within the framework of 



the relativistic quark model based on the quasipotential approach. This is different from the 
discussions based on the quasipotential approach with the same structure of heavy quark 
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mass corrections predicted in the usual HQET |L7j, where the decay rate of B — > D^lP is 
only slightly increased by subleading 1/mg corrections. 

In comparison with the experimental data reported by CLEO [|18|] and ALEPH [19| 
groups, our result for the branching ratio of the B — > D\lv decay with the inclusion of 1/mg 
corrections is in agreement with both measurements. While for the B — > D^lu decay, the 
result at itlq — ► oo limit is within the CLEO upper limit but not within the ALEPH one. 
When including 1/mg contributions, we find that the resulting decay rate for Br(B — > D^lu) 
seems to go over the CLEO upper limit though it may still be consistent within the large 
erros due to the big uncertainties of the choices of the thresholds. On the other hand, in 
deriving the results in Table 1, we have made the assumption that rff = 0. In general, 
the contributions of rff may not be neglected. Their effects have shown in Fig. 4, where 
we have used the above sum rule results for r, t/q and r/g and the thresholds u c = 1.9GeV, 
u** = 2.95GeV, u* = 2.35GeV and = 1.85GeV. We have also made the assumption 
that rff and f have the similar dependence on y, f « f (l)/(l + ^i) with a 2 « 0.7. It is seen 
from Fig. 4 that rjf, 7/2, 7/3 and rj b influence the branching ratios in different ways. When r]\ 
becomes larger, BriB — > DiW) decreases while BriB — > B\lv) increases, and 7/2 influences 
the two branching ratios in the same manner as r\\ does, but both branching ratios are not 
sensitive to r^. 7/3 has opposite effects on the two decay modes, i.e., its increment enlarges 
Br(B — > Dili?) but suppresses BriB — > Dl^v). For r/ b , it always appears in the form 
so it affects Br(B — > D\lv) and BriB — > D^W) in the same way, i.e., a negative value of 
if may suppress both branching ratios. Some reasonable results for the branching ratios at 
certain values of rff are listed in Table 2 and Table 3. From Table 2, Table 3 and Fig.4, 
we see that both branching ratios of B — > D^D^lv may be suppressed after considering 
the possible effects due to the contributions of chromomagnetic terms at 1/rriQ order. The 
resulting two branching ratios can easily be made to be consistent with the experimental 
measurements when r/ 6 , r/J, r]\ and r/3 are in the proper ranges. 



V. SUMMARY 

The HQEFT has been reviewed and applied to study the semileptonic B decays into 
excited charmed mesons (Di^D^) with jf = | + . The form factors of the matrix elements 
relevant to these two decays have been expressed in terms of wave functions within the 
framework of HQEFT. It has been shown that with the inclusion of quark-antiquark coupled 
sectors, the relevant matrix elements can be parameterized by the leading order Isgur-Wise 
function and additional twelve wave functions rff [Q = b,c; i = 0, 1, 2, 3) and Kj (j = 1, 2) 
appearing at l/mq order. Two wave functions rf and rf characterizing the l/rriQ corrections 
of effective current in the usual HQET |I],|2| are absent in HQEFT. 

By adopting proper interpolating currents for the heavy-light mesons, the leading order 
wave function r and two important wave functions i]q and t]q at 1/rriQ order have been 
calculated by the QCD sum rule approach. Zero recoil values of k[ and k' 2 , which are 
1/mg order wave functions for the matrix elements between | + excited states, have also 
been extracted from fitting the excited meson masses. By using these calculated results 
and also considering possible effects from wave functions rff (i = 1,2,3) arising from the 
chromomagnetic operators at 1/rriQ order, B — > DilD^lv decay rates and branching ratios 
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have been evaluated. When neglecting rff (i = 1,2,3), we have shown that both of these 
decays receive large 1/rriQ corrections from t]q and k\ (1) (i = 1,2). For B — > D\lv decay 
this is similar to the results obtained from quasiquark potential approach based on the same 
structure as in the usual HQET, but for B — > D* 2 lv decay the situation is quite different from 
the case described by that approach, where BriB — ► D^lv) only increases slightly when the 
1 / tuq order contributions from are included. It has been shown that when considering 
all possible contributions from wave functions at the I/tuq order in the HQEFT, but without 
breaking the I/tuq expansion, the resulting branching ratios of the semileptonic B decays 
into | + excited charmed mesons can agree well with the experimental measurements for the 
proper ranges of rjf . 
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APPENDIX A: LOCALITY OF HQEFT 



For completeness and clarification, we briefly review the formulation of heavy quark effec- 
tive field Lagrangian that keeps both effective quark and antiquark fields ||. In particular, 
we would like to emphasize, by adding this appendix, that the HQEFT is a local effective 
theory and contains no non-local operators. 

Firstly, denote the heavy quark (antiquark) field as 



Q = Q+ + Q-, 

with Q + and Q~ corresponding to the two solutions of the Dirac equation. 
Defining 



(Al) 



Qt 



- - ^ 



R„ — 



(A2) 



with v M an arbitrary four- vector satisfying v 2 = 1. Furthermore are related to the desired 
effective fields defined as 

Qv 



Qv: Qv Qv^ 



—ifrriQV-x 



(A3) 



From the equation of motion of quark field and antiquark field, the original fields Q and Q 
can be expressed by the new fields as follows, 



G ± = [l + (l-^)- 1 |fe]* : =^ 

= e T* mfl - [ l + (l_^ ) -l^]Q± = 

Q ± = ^[i + ^(i - 4£^n = Q v * 



(A4) 



Qti± + 



2m,Q 



Q 



:i- 



2rriQ 



-lj e ±im,QV-x = Q± e ±irriQV-x 
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D^, Ip\\ and Ip± are defined as 



JkD' 1 <p=-J kD^, 
k = jp-j{y-D). 



The QCD Lagrangian becomes 



QCD 



flight + ^ 



where Cu„ht represents the part of Lagrangian containing no heavy quarks, and 



C 



eff 



r- 



C: 



c; 



with 



zmn zmn 



2m 



Q 



Am 2 Q 



2m Q 



2m Q 

x (i - %±^r' iW « . &m 



2rn Q 



-*#0(i 



2m Q 



x (1 - Jij-I^JQJ = e ±2 ^^Q±i?gj 
2m Q 



(A5) 



(A6) 



(A7) 



(A8) 



It is seen that all parts of the effective Lagrangian are local. 

When quark fields and antiquark fields decouple completely, it is reasonable to deal with 
only section £+ + or L~~ independently. This is just the case considered in the framework 
of the usual HQET. To consider the finite quark mass corrections one should also include 
the contributions from Cf T and L~~ as well. 

Similarly, the heavy quark currents can also be decomposed into four parts in a similar 
way, 

J(x) = Q'{x)TQ{x) = Q' + TQ + + Q I+ TQ- + Q'-TQ + + Q'TQ- 

j eff (x) = gy S TuQ+ + 4',+ S tqq- + q'j £ r^g+ + gy S r^g; 



7 +H 



JT + ^ + JT- (A9) 

The matrix element on the r.h.s. of fl2.2| ) is actually evaluated via 
(H'J J eff e^ c ^\H v ) = (H' V ,\(J+ + + J+- + J~ + + J;-)e l ^ ct++ct ~ +c ^ ++c ^\H v ). (A10) 
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Here J e ff and £ e ff shall include all 4 parts instead of only the ' + +' parts in ( |A7|) and ( |A9|) 
as has been done in the usual HQET. The first sector ' + +' of J e ff contributes 

(H'JQ^ S TLuQte l f c *"\H v ). (All) 
If neglecting all except ' + +' sectors in J e ff and £ e //, we get 

(H'JQp&ruQ+e'lftXttlHv), 

which is just what the usual HQET treated. 

The contributions from the ' H — ' sector of J e // is 

(H'JQp S TG}Q-e^ c 'ft\H v ). (A12) 

In the case that both the initial and final states contain only heavy quarks (no heavy 
antiquarks), matrix elements such as 

(H' vl \^ZTuQ-\H v ) 



do not contribute, and (|A12 ) contributes only at higher perturbation order, i.e., only when 
C e ff is inserted into the matrix elements. Therefore the leading order contribution from 
( gig) should be 



(H' vf \i J dxdyT{(Qp I TtuQ-)(x), (Q-BQ+){y)}\H v ) (A13) 



and then the effective antiquark fields Q v (x) and Q v (y) should be contracted. Due to (|A8|) . 
this will yield the propagator 



-p + m Q + 0(l/m Q ) 



which is of 0{l/rriQ). 

Contributions from other sectors of the current J e f f can be treated in the same way. It is 
easy to see that (H' v ,\ J+ + e l S Ceff \H v ) gives contributions of 0(1), and (H' v/ \J± T e l S Ceff \H v ) 
is 0(l/mg), whereas (H' v ,\ J~ v e' 1 -f Ceff \H v ) is 0{l/mV) since C\/ mQ should be inserted twice, 
and each contraction of Q~ and Q~ gives a 1/toq suppression. 

To be more clear, contracting Q~ and Q~ in (|A13|) yields 



(H'Ji / d^yQ' v f{x) Q Too / -\- ~ , nh , , BQ+(y)e^ C ^\H V ). (A14) 

J J {27ipfv -p + m Q + 0(l/m Q ) 

which can be written in the following form with replacing the momentum p by the derivative 
d and performing the integral of the momentum 

(H'Ji [ d 4 xd 4 y^(x) I TS(x - y)u - \ -BQ+(y)e l I C ^\H V ), (A15) 

J ifv ■ o + mq + (J{l/mQ) 
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Using the same trick, one can treat all contributions from (H' V ,\J+ e l f Ceff \H V ) and obtains 
{H' v ,\J^-e^ c ^\H v ) 

= (H'JQ'J tJ Tu[-i% - ^-0±O- - ^r l J ]- 1 BCfte i S c «f\H v ) 

= (H' v ,\e^ m Q' v '- m ^ v > x Q^ tj Yu{-A- x )BQ+e^ ~ Ceff \H v ) (A16) 

which means that, effectively, one can reexpress J + ~ to be the effective current in which 
only the effective quark fields Q'£ and are used. Namely, 



J+- = (Q + TQ-) = to VujQI 

- $; + S VCj(-A- 1 )BQ + u = j(™' Q ,v'-™QV>*Q>+ Zj Tu{-A- l )BQ+ (A17) 



In an analogous way, one can reexpress J„ + , J v , £t~-> £-v + an d £-v~ the corre- 
sponding effective currents and Lagrangians by only using the effective quark fields Q'$ and 
Qv- Consequently, we have 

£ ~^ £eff — ► £ e ff = £t + + £t + -> ( A18 ) 

2m Q 2m Q 2m Q z^)|| 

x [1 - ^ + ^-iy^m^^rWl + ^-^±(1 " ^L)- 1 ^^ 
2m Q 2m Q %Jp\\ %JP\\ 2m Q 2m Q 

1 -Q+A^-iip^l - i^L + J_^ J __L^ J _]-i^_ L _l-AQ+ (A19) 



2mg " 2mg 2mg i^y 

which represends the additional contributions to the effective Lagrangian £+ + that has been 
widely adopted in the usual HQET. This additional part may be regarded as the effective 
potential part of heavy quark due to the exchanges of virtual antiquarks. It is seen that 
when one imposes the on-shell condition AQ+ = 0, i.e. £+ + = 0, the effective potential part 
£+ + also vanishes, i.e. £+ + = 0. For off-shell case ~ (Hp±) ~ A, or 

~2^ <<c w> 

the leading contribution of the effective potential part is 

Correspondingly, the heavy quark current turns out to be 

J = Q'TQ — > J eff — > J+ f + f = J+ + + J+ + 

= e^'-^-^Q'+TQ: + J-Q'+T^-(zip ± ) 2 Qt 

2m Q iJP\\ 

2m Q' -i m Am Q ww 
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i 



Am Q >m Q 



g;t(-^ ± ) 2 x 



1 -r4r(^) 2 ^ + 0(^)} 



#11 



^11 



(A20) 



with 4% the leading term jf} f = e i{ ^Q' v '^ v > x Q^YQ+ and J^/ f mQ) the remaining terms 



in J, 



eff- 



APPENDIX B: LORENTZ STRUCTURES 

Here we present the general Lorentz structures of Ly^ a ^ appeared in eqs. (p. I|)-(p76| 
1 



VE 



(y 2 - \)f v + 3v^v" + (1 - 2y)vV - 2>yvN v - (y - 2)v'»v v 



£% = i—(l+yy™Pv a v' /3 , 



c 



fia/3 
2V 



_i (e^P{ Vp - v' p y) + e^(v a - v' a y)) v a v' p 



(Bl) 
(B2) 
(B3) 



Of = \(i + y)g aP (v» - v'n - \{l + y)g a ^ - xff>y) - 1(1 + y)g^(v a - v' a y) 

+ -(1 - y)v' a v p v^ + -(1 - j/^W" - -(1 + W 
2 2 3 

- ~(1 + y)v' a v' p v'^ + (v'V 3 + uV)t/". (B4) 
o 
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Fig. 3c 



Fig.3. Variations of the wave functions r, t/q and t]q with respect to y at T = lGeV. 
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t,-,(1) (GeV) 

Fig.4c 

Fig.4. Branching ratios of B — > D^D^) decays at different values of rjf . a. r/ b (l) = — 0.4GeV 2 ; 
b. n b (l) = -0.6GeV 2 ; c. r/ 6 (l) = -0.8GeV 2 . In these figures 77^(1) = -0.6GeV 2 is used, and we 
assumed that rff and f have the same monopole. 
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Table 1. Decay rates F (in units of |Vd,/0.038| 2 x 1CP 15 GeV) and branching ratios BR (in 
%) for B — ► Di(D2)lv decays in the infinitely heavy quark mass limit and taking account of first 
order 1/ron corrections. R is a ratio of branching ratios taking account of 1/mg corrections to 
branching ratios in the infinitely heavy quark mass limit. 







5 -» Di/P 


B -» D^z? 


rriQ — >■ oo 


r 


1.87 ±0.53 


3.24 ±0.94 




Br 


0.45 ±0.13 


0.79 ±0.23 


with t]q 


r 


2.42 ±0.65 


4.05 ± 1.09 


(') 
and k\ 


Br 


0.59 ±0.16 


0.99 ±0.26 


R 


1.30 ±0.54 


1.25 ±0.05 


Experiment 


Br(CLEO) H 


0.56 ±0.13 ±0.08 ±0.04 


(0.8 




Br(ALEPH) j|| 


0.74 ±0.16 


(0.2 



Table 2. Branching ratios (in %) at some values of rjf. 772(1) = — 0.6GeV 2 and the monopole 
assumption are used in deriving the branching ratios in this table. 



r/ fe (l)(GeV 2 ) 


r/f(l)(GeV 2 ) 


^(l)(GeV 2 ) 


Br(B -> D x lv) 


Br(B -» D%to) 


-0.4 


-0.2 


-0.1 


0.70 


0.52 


-0.2 


0.1 


0.80 


0.41 





0.2 


0.47 


0.58 





0.4 


0.58 


0.46 


0.2 


0.4 


0.32 


0.71 


0.2 


0.6 


0.42 


0.58 


-0.6 


-0.2 


-0.1 


0.66 


0.45 


-0.2 


0.1 


0.75 


0.35 





0.2 


0.43 


0.51 





0.4 


0.53 


0.40 


0.2 


0.4 


0.29 


0.63 


0.2 


0.6 


0.39 


0.51 


-0.8 


-0.2 


-0.1 


0.61 


0.40 


-0.2 


0.1 


0.71 


0.30 





0.2 


0.40 


0.45 





0.4 


0.49 


0.35 


0.2 


0.4 


0.27 


0.56 


0.2 


0.6 


0.35 


0.45 
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Table 3. Branching ratios (in %) at some values of r}f . 773(1) = OGeV 2 and the monopole 
assumption are used in deriving the branching ratios in this table. 



r/ 6 (l)(GeV 2 ) 


) 


77 3 (lJ(GeV ) 


oryD — > U\iU) 


o„/ d , D*lrA 

DVyD — > D^W) 


-0.4 


n 
— U.z 


n 1 
— U.l 


U.00 


U.Ol 


n 
—U.z 


n 1 
U.l 


U. (o 


u.4y 





U.z 


n /1 
U.4z 


U.b 7 





n A 
0.4 


U.OZ 


U.00 


n 
U.z 


U.4 


n on 
U.zy 


U.ol 


n 
u.z 


n p. 

U.D 


U.00 


U.D ( 


-0.6 


—U.z 


n 1 
—U.l 


U.Oo 


U.04 


n 
—U.z 


U.l 


U.bo 


U.4o 


U 


U.z 


U.00 


U.OU 


n 



U.4 


U.48 


n A q 
U.48 


0.2 


0.4 


0.26 


0.73 


0.2 


0.6 


0.35 


0.60 


-0.8 


-0.2 


-0.1 


0.54 


0.48 


-0.2 


0.1 


0.63 


0.37 





0.2 


0.35 


0.53 





0.4 


0.44 


0.42 


0.2 


0.4 


0.24 


0.66 


0.2 


0.6 


0.32 


0.53 
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